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Abstract 

The BK equation in the conformal basis is considered and analyzed. It is shown that at high 
energy a factorization of the coordinate and rapidity dependence should hold. This allows to simplify 
ON ' significantly the from of the equation under discussion. An analytical ansatz for the solution to the 
BK equation at high energies is proposed and analyzed. This analytical ansatz satisfies the initial 
condition at low energy and does not depend on both rapidity and the initial condition in the high 
energy limit. The case of the final rapidity being not too large is discussed and the properties of the 
transition region between small and large final rapidities have been studied. 



X 

- 1 Introduction 



The scattering of two distinct object, such as hadron and nucleus, or DIS at high energy is described 
by so-called " fan" diagrams where only splitting of one pomeron into two is taken into account [H GJ [3] . 
This situation was already discussed many years ago in the framework of the phenomenological pomeron 
[1], and some early attempts have been taken to generalize it for the case of QCD [5]. The equation 
for the "fan" diagrams in the operator expansion formalism was written by I.Balitsky [9], and in the 
dipole model framework [11] by Yu.Kovchegov [TO]. The resulting Balitsky-Kovchegov (BK) equation is 
very well numerically studied in both, the saturation region and the region of small non-linearity (see 
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[T2| [T3[ [H] and [J5J HE])- One of the important features of the BK equation is the presence of the 
saturation scale at high energies that increases exponentially with rapidity, and a geometrical scaling of 
the solution [TTl [TBI HH] ■ Despite the well understanding of the properties of the BK equation, the question 
of an exact analytical solution to the BK equation still remains open. This full analytical solution may 
be very useful in the applications of the BK equation to different scattering processes at high energies. 

In the present paper we want to bridge the gap in the analytical study of the BK equation. We 
consider the formulation of the BK equation in the conformal basis and solve the problem using some 
simplifications. In the conformal BK equation we keep only the terms with zero conformal spin, which 
corresponds to leaving only the BFKL pomerons propagators in the "fan" diagrams. Due to highly 
complicated forms of the integrals appearing in the calculations the exact analytical solution is beyond 
the scope of the present paper. However, we do discuss the possible ways of finding this exact analytical 
solution. We propose an ansatz of the solution to the BK equation, which relates between high energy 
behavior and the given initial condition at zero rapidity for the pomeron field. This ansatz we call 
" phenomenological" or the BFKL ansatz for the solution to the BK equation. The reason for this is very 
simple, "phenomenological" because this solution is similar to phenomenological fan solution found in [4] 
(see also [20J for more details) and the BFKL since the found solution corresponds to the leading BFKL 
terms in the "fan" structure leading at high energies. 

The said ansatz was obtained using the conformal invariance of the pomeron field and the action of 
the theory at high energy. The conformal invariance disappears when the final rapidity becomes not too 
large. This happens due to the influence of the non-invariant source that cannot be neglected at low 
energy. Therefore, for small final rapidities we consider a similar ansatz, which nevertheless does not 
possess property of the conformal invariance. Having two solutions for two different regions we are able 
to analyze the properties of the transition of the pomeron field from one region to another. 

The present paper is organized as follows. In the next section we consider the BK equation in the 
usual coordinate formulation and rewrite it in the conformal basis. In Sec. Owe consider the conformal 
invariance of the theory at high energy and obtain scaling properties of the pomeron field. In Sec. 0] 
we propose a possible ansatz of solution to the BK equation in the conformal basis at high energy and 
discuss the analytical structure of the solution. In Sec. [5] we consider an ansatz for the pomeron field at 
small values of rapidity. Sec. [6] is dedicated to the conditions and scales arising in the region where the 
transformation of the pomeron field from the region of low energy to the region of the high energy place. 
In the last section the conclusions and discussions are presented. 

2 BK equation in conformal basis 

We begin our consideration with the effective field theory of the interacting pomerons. It was shown by 
Braun ( see Refs. |2H [22| [25} 124] ) that one can describe the interacting pomerons in the large N c limit 
in terms of the effective action. The effective action with pomeron splitting only reads 

S = Sq + Si — Se , (1) 

where Sq is a free part of the action 

So = J dydy' d 2 p 1 d 2 p2d 2 p[d 2 p 2 ^(y,p 1 ,p2)G y \{p 1 ,p2\p[,p 2 )^(y,Pi,P2) ( 2 ) 
Si is a interacting part of the action 

a 2a 2 N c f f d 2 pid 2 p 2 d 2 p 3 t , 

Si = — / dy / K 2 2 2 ( L i3 $ T (y> Pu Pa)) *(y, Pi, Pa) $(y, P2, Ps) (3) 

n J J r 12 r 23 r 31 
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and Se is a source term of the action 

Se = J dy J d 2 p 1 d 2 p 2 <S>(y,p 1 ,p 2 )T A (y,p 1 ,p 2 ) + &(y, Pi, P2) tb(v, Pi, P2) (4) 
The propagator of the theory is defined through the BFKL Hamiltonian [27\ [28j [29] as follows 

G-\,(pi,P2\Pi,P2) = (v^V? (A + # (pi)P2 )^ s 2 ( Pl - p' 1 )5 2 (p 2 - p' 2 )d(y-y) (5) 

The source terms in Eq. (JH) are the initial forms of the functions $ and <&t at rapidities y = and y = Y 
respectively 

${y,pi,P2) y =o = T B (y,Pl,P2) = f B (pi, P2) S(y) (6) 

&(y,Pi,P2)y=Y = T A (y,p 1 ,p 2 ) = f A {pi,p 2 )5(y - Y) (7) 

The pomeron field <£(y, ri,rj) in Eq. ([2]) and Eq. Q is the generalized (skewed) gluon (parton) distribution 
written in coordinate representation, see |25U26| . and L13 is the Casimir operator of the conformal group 



Lis = pWivl = pf 3 V?V| (8) 
In this formalism the BK equation is the equation of motion for the pomeron field &(y, Pi, pj) 

(JL + Hfapa)) $(y,pi, P3 ) + / ^4 L Hy,pi, P 2)Hy,P2,P3) = (9) 

with the initial condition for the field given by Eq. ([6]). 

For our further analysis it is more convenient to write Eq. Q in the conformal basis of functions 
E fM (p 1 ,p 2 ) (see Ref.[27]) 

E,( Pl ,p 2 ) = ( -^-) h ( h (10) 



P10 P20 / \ P10 P20 

where p = {h, po}> h = -Mp + iv and h = 1 — h*. For the sake of simplicity we adopted notation used in 
Ref . |23j . These functions are the eigenfunctions of L13 operator 

Lv&Epfapz) = \- l E^ Pl , P2 ) (11) 

with the eigenvalues 



1 

[{n + l) 2 + Au 2 J[(n - l) 2 + 4v 2 } 
The pomeron field can be expanded in this basis as follows (see [27] 



®(y,pi,P2) = E Mpi>P2)Mv) ( 13 ) 

where 

My) = I p \ P2 E;( P1 ,p 2 )$(y, P1 ,p 2 ) (14) 
j P12 

In Eq. ()13p and hereafter by the conformal summation one should understand 



H n=—oo 



E = E U**-^}*!* ^ 
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Recasting the Lagrangian of the theory in the conformal basis ( for more details see [26J) one obtains 



dy 



- U^y) *t (y) - V A , ttlV $+(y)^(y)$,(y) J + r MA + *t (y) r MB | (16) 

with u;^ being the BFKL eigenvalues 

/ I n I + 1 \ 
o;, = W/l = aU(l)-Re^ l - +iv)\ (17) 

and withe the sources Tu A and r^ B for the fields and $ at rapidities y = Y and y = correspondingly. 
The effective action in this form leads to the equation of motion for each field 

= »**M - ^ E *m(v)*/«(v)^i.i^ ( 18 ) 

where V^ jAtljAt2 is the triple pomeron vertex in the conformal basis given by 

d 2 pi d 2 p 2 d 2 p 3 
Pxi P23 P31 

and jj, = {1 — /i, po}. The expression of Eq. (fT8l) is, essentially, the BK equation in the conformal basis 
and its solution is the main subject under consideration in the present paper. 

Further simplifications of Eq. (|19p using properties of the conformal group were performed in Ref.[31j 
and the simplified expression reads 

V wm =n(h ,h 1 ,h 2 )(z -z 1 )- Aoi (zi-z 2 )- Al2 (z2-z 3 )- A ™(zo-z 1 )- A n (20) 

where A12 = hi + h 2 — ho and A 12 = hi + h 2 — ho. We follow the notation of Ref.[30l [3Tj and denote 
po of po by Zo, i.e. po = {ho,zo} in Eq. ([20]) and below. Using the simplified expression for the triple 
pomeron vertex in Eq. (|20p we can write explicitly the nonlinear term in Eq. (|18p as 



Vwuv, = I d P 2 d 2 2 d 2 P3 E^px, pa) (p 2 , Pi) E m (p 3 , Pl ) (19) 

J Pi 9 Am P.^i 



1 9™ 2 AT 00 /*oo „2 i 5i r oo ,,2 _i_ !il /• 

" " ' E / ^1^^ / -^T^ / ^1 / (Zo ~ Zl )-^n {Z1 - , 2 )-A 12 (, 2 _ Zo) -A 2fi (21) 

, J — 00 TT J — 00 TT J J 

x dzi dz 2 (z - zi)" A51 (2i-z 2 )~ Al2 (^2-2o)" A2 o$ Ml (y,zi)^ 2 (y,22)0(l-/io,/ii,/i2) 
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ni ,ri2=— 00 



where the A^ notation stands for /io to be replaced by 1 — ho in the expression for A01 defined above. 
This comes from the fact that the expression for Eq. (|19p given by Eq. (|20p appears in Eq. (|18p with 
one of the conformal vertex functions conjugated. To see the scaling properties of this nonlinear term we 
introduce dimensionless variables 

wi = — , w 2 = — (22) 

zo z 

and their conjugate. In terms of these new dimensionless variables Eq. (|21|) reads 



2ajN c 

TT 



E 



ni,n 2 = 



dv\ 



TT'' 



di>' 



14 + 



TT'' 



z 



*0 
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x J dwi J dw 2 (1 - u;i)~ A ei (w 1 -w 2 )~ Al2 (w 2 -l)~ A ^o J dm J dw 2 (l-wi)~ Adl (m - w 2 )~ Al2 (w 2 - l)" Aa 

x $^ 1 (y,z wi,z Q w 1 )®^ 2 (y,z w 2 ,z w 2 )tt(l -h ,hi,h 2 ) (23) 

As one can see from Eq. (I23p the scaling properties of the triple pomeron term are merely determined 
form the properties of the pomeron field and the power of the dimensionful variable zq. 

3 The scaling property of the pomeron field 

In this section we show how the scaling properties of the pomeron field suggest an ansatz for solution 
of the conformal BK equation at high energies. We come back to the scaling independent value of the 
theory, to the action given by Eq. (|16|) . and check the scaling properties of the pomeron field. We assume 
that our theory is conformal invariant at high energies. As a consequence of this fact the pomeron field 
scales accordingly to its dimensions 

resulting in the following scaling property 



[%(Y, Zo )] = zt 1 zt 1 (24) 



%(Y,Xzo,Xzq) = A A ®n(Y, zo,z ) (25) 

Next we analyze the scaling property of the source. We plug the equation of motion Eq. (|18[) into the 
action of the theory and obtain 
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dy ®M r, A (y) = E / dv / d " *> W «>) f ^o) , (26) 

fx n=-oo J -°° n J 

where Y is the final rapidity of the process and f flA (zQ) is defined by 

Wy) = fpMSiy-Y) (27) 

Using the scaling property of the pomeron field found in Eq. (|25|) one can derive the scaling of the 
source 

\ h \ h f^ A (\z ,\z ) = fp A (zo,zo) (28) 

Another consequence of the scaling property of the pomeron field Eq. (|25p is that it preserves unchanged 
the form of the equation of motion with the triple pomeron term given by Eq. (12311 . This observation 
and the condition of the invariance of the action under the rescaling prompts the following ansatz for the 
pomeron field at high energy 

^ lx {y,z) = z h - 1 z }l - l f^\y) (29) 

It is easy to see that this ansatz satisfies all conditions which we have on the theory at high energy. An- 
other important observation we made is the fact that Eq. (|29j) presents a manifestation of the holomorphic 
separability of the pomeron field at high energies. 

Now we are in position to plug the ansatz Eq. (|29p into Eq. (JTHJ) with the help of Eq. (|2ip and read 
out the relevant integrals over the coordinate variables. The coordinate integration factorizes and is 
performed in the Appendix A. The resulting equation is the integro-differential equation for f^ h ' h \y) 
function given by 

df{h d h)(y) = Wft /("■%) E f (hl '~ hl) (y)f ih2 ~ h2 \yMi- Mi, h 2 )i{\-hMM)i£-hMM) (30) 

hi,h 2 
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where the integrals 1(1— h, hi, h 2 ) and 1(1— h, hi, h 2 ) are calculated in the Appendix A and the summation 
over h denotes 

The solution of Eq. (|30f) would give an exact analytical solution to the BK equation in the limit of high 
energy. 

It is interesting to note that the form of the coordinate dependence of the function <3?^ (y) given by the 
ansatz Eq. f|29[> resembles the form of the momentum dependence of the solution of the BFKL equation 
for zero transferred momentum. The physical meaning of this dependence is different since the variables 
z and z are conjugate to the transferred momentum while in the solution of the BFKL equation such 
power dependence arises in the transverse momentum of the pomeron. However, the question of a more 
subtle relation between the two cases is still to be answered. Another interesting observation, is that a 
factorized form of the ansatz Eq. ()29[) is a consequence of the fact that at high energy the solution of the 
equation of motion does not "remember" about the initial conditions and there is no another scale which 
can destroy this factorization. 



4 Solution of BK equation in conformal basis at large final rapidities 



As it was already mentioned Eq. (I30p is our master one variable integro-differential equation, solution 
to which gives the full analytical solution to the BK equation at high energy. Unfortunately, being 
much simpler than the original BK equation it still quite complicated. The exact solution of Eq. (|30l) 
is beyond the scope of the present paper and we only wish to find the leading solution at high energies. 
We consider the case of zero conformal spin ni = which corresponds to the fan diagrams with only the 
BFKL pomeron propagators. For this case Eq. (|30f) reads 

df v {y) _ fVI , a 2 s N c '<*> - 2 '°° - 2 



■^VT J — 00 J— oo 7T 

n f 1 • 1 • 1 • \ r2, 1 .1.1., 



dy 



We want to investigate the properties of the solution to Eq. (|32|) assuming that the function f u (y) can 
be expanded in powers of v 

00 

r(y) = J2My)» k ( 33 ) 

k=0 

and plug into Eq. (|32j) . To the zeroth order we obtain 

9 My) . ( v a 2 s N c 2 /- u 2 /- u 2 

11 1 . \ .,1 1 1 ...,11 1 



x ^ ( 2 + ivu 2 + ^2' 2 + 2 + iv2)I ^ 2 +iUl, 2 + 



or the same equation in more compact form 



dy = ^ojo(y) — ^— j {y)C (35) 



with some coefficient 



00 v 2 j-oo y 2 (11 1 \ 1 1 1 

C =l ZJ^I ■^du 2 Vt[-,-+iv 1 ,- + iu 2 \l 2 (-,-+iui,-+iu 2 ) (36) 



-4 — L / ~4 

00 " J —00 " 
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The solution to Eq. ([35]) is readily found 



My) = n 2 N n ( 37 ) 

where the unknown coefficient Coeff is to be determined from the initial conditions for fo(y) function. 
Let's now simply write this condition in the following form 

f (y = 0) = F in (38) 

which leads to 

, / \ exp(u> y)F in 

My) = -tpn F~, — ; — : : • ( 39 ) 

^F in ^ (exp( Wo y)-l) + 1 

We see, that at asymptotically high energy this ansatz does not depend on the form of F{ n , leaving the 
possibility for arbitrary form of F{ n . We will precisely define this function latter. For obvious reasons 
we can call Eq. (|39p the double leading ansatz to the solution of the BK equation. First, because we 
consider only the leading term dominant at high energies, which is the BFKL fan structure; second, the 
function fo(y) is the first term in the expansion in the powers of v . For the reasons we can call Eq. (j39H 
the BFKL ansatz for BK equation. This ansatz is also of a great interest since it is similar to the solution 
to the analog of the BK equation in the phenomenological pomeron theory, see for example [20J. 



As a next step in our analysis we partially reconstruct v dependence of the function f u {y). The form 
of Eq. (f3U|) suggests the following form of the solution with explicit v dependence 

rw= F ... (40) 

Fin % (y) ( exp(w A1 y) - 1 ) + 1 
where go(is) is the first term in the expansion of f(y) in inverse powers of one pomeron exchange 

oo 

r(y) = Y,to(y) e ~ kUvV - ( 41 ) 

k=0 

The function g n (u) are found from the equation Eq. (I32p in the form of 

°° ^2 at poo ,.2 poo ,,2 

/i i i \ 2/ i i l 

As it was already mentioned due to complexity of v dependent functions in the kernel of this equation 
we are only able to find the solution leading in energy. This means that we keep only terms with k, k\ 
and k 2 equal to zero obtaining equation for go (v) only 

n 2 AT roo ,,2 poo ,,2 

^!»W = -^ du 2 ^ du 3 ^g (u 2 )g (u 3 )ni 2 (43) 

J — oo 7T J — oo 7T 

Plugging Eq. (|40|) into Eq. (|29|) we obtain, that when final rapidity of the process is large, the asymptot- 
ically leading solution of BK equation may be represented with the help of our ansatz 

•„M = M- 1+2 *" _ ^f F '~ (44) 
Fin g iy) { exp(u^ y) - 1 ) + 1 



7 



As it must be, at the asymptotically large rapidity the pomeron field Eq. (pHj) can be written as 



=\z \- 1+2 ^g (u), (45) 

that confirmed the scaling property of the pomeron field at large rapidities. In the case of the rapidity 
smaller then the final rapidity Y of the process, the factorization determined by Eq. (|29p is already 
not necessary and may be broken. The possibility of the scenario when the final rapidity Y is small is 
discussed below. 

The ansatz Eq. (|4U|) correctly reproduce the property of BK equation at large rapidities, namely the 
independence of the ansatz on the initial conditions, the constant behavior at high energy and scaling 
independence of the action for the pomeron field. Nevertheless, the function Eq. (|40p being used in the 
amplitude only at rapidity y = Y, could be defined at all values of rapidity. We find the form of _Fj n 
function from Eq. (|44p taking y = in the ansatz 

^(y)y=o =\z \- 1+2iu F m (46) 

As it was mentioned in the first section of the paper, the Lagrangian of the theory includes the source 
terms. This implies 

Mv)8(y) = Tp B 6(y) (47) 

resulting in 

\zo\~ 1+2lu F in = f UB (48) 
Plugging Fi n from Eq. (I48|) back into Eq. (I40p we obtain for our ansatz 

M 1 * %v Tp B g {v) ( exp(tJ M y) - 1 ) + 1 

The analytical ansatz of Eq. (149)) interpolates between two desirable features of the behavior of the 
solution to BK equation. At high rapidity it is determined by the function which does not dependent 
on rapidity and initial condition of the problem, whereas at rapidity zero it is equal to the given initial 
function of the pomeron field. It is important to underline therefore, that the source field f UB in Eq. (|49p 
is arbitrary due these properties of the ansatz. It is washed out at high rapidity and there are no special 
constraints on the functional form of this source. This source function may be arbitrary and may depend 
on some external scales of the problem. 



5 Solution to BK equation in conformal basis at small final rapidities 

In this section we consider a solution at small final rapidity Y. In this case the asymptotic expansion 
in powers of e~ Utiy in Eq. (|44p is not valid anymore. The proposed ansatz given in Eq. (]49p does not 
describe correctly the solution to the BK equation and does not provide the scaling invariant solution for 
the amplitude as well. It is not surprising since we do not expect that at low final rapidity the solution 
will preserve this scaling invariance property. Indeed, let us return to the triple pomeron term in the 
equation of motion Eq. (|23|) . keeping only zero conformal spins in the formulae. Clearly, all our previous 
consideration are valid if we can justify the asymptotic expansion of the pomeron field in the form 

oo 

Mv) = Es*W e " Hs ( 5 °) 

fc=0 

However, if Y is small another expansion should hold 

oo 

MV) = (51) 

k=l 
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Similar the previous case we put the expansion of Eq. (|5ip into the equation of motion keeping only the 
first term of this expansion. It is easy to see that at small rapidities the triple pomeron term is not 
enhanced by rapidity exponential and can be safely neglected. In this case we obtain the BFKL pomeron 
solution which corresponds to the first term in the expansion Eq. (15 1|) 



®»(y)=f m e^y (52) 

and determines the form of the action 



TfM B (Xz ,\z ) = A 1+2 w t^ b (z ,z ) (54) 



/OO j,2 r 

dv— / d 2 z Q f IXB {z Q )e w ^f llA {z Q ) (53) 
-oo 7T J 

Similar to the case of the pomeron field we deduce from the from of the action the scaling property of 
the f Ms , namely, 

'' p-b \" ^uj " ^u; — <t >jj, B \ 

Naturally, for an arbitrary source field f^ B such a strong constraint, in general, is not satisfied. Therefore, 
we do not expect the invariance of the action under the rescaling of the variable zq at small energies. 
Such an invariance restores at large energies, when the dependence of the amplitude on the source f^ B 
disappears. 

As a next step in our discussion we consider the first two terms in the expansion Eq. (|5ip of the 
pomeron field at small energy 



My) = (n^o) - g?(zo))<r»y + - g ^)e 2 ^ y . (55) 

Using equation of motion we obtain the equation for the <Ji(zo) function from this expansion 

#fo)e^» WM = -^?£ £ V~ m2 (f^ B { Zl ) - fff(*i))(f M2fl (* 2 ) - 9t 2 (z 2 ))e^ + ^y (56) 

In this way we can consider an expansion with any arbitrary number of terms obtaining a chain of 
the equations similar to Eq. (|56p with involved interference terms between the functions g^(zo). Instead, 
we find an ansatz that have a structure of the expansion Eq. (|5ip and will coincide with the expansion 
Eq. (151ft to some order. As the simplest example we take expansion up to the second order and the form 
of this ansatz is borrowed from Eq. (|49|) as follows 

<j> ( v) = expK^f^ f57) 

" W l^i-si" r Ms 50(^0) (exp^yl-l) + 1 1 ' 

Expanding Eq. (I57p and comparing the first two terms of this expansion to the Eq. (|55j) we obtain 

g^z ) = -\z \ l - 2 ^f 2 B g (v,z ) (58) 
We use this expression in Eq. (|56p to obtain the equation for the function go(v, Zq) , namely, 

\zo\ l - 2iv fl B go(v,z ) e 2 ^u, M = 

E y^^ B {zi)f, 2B {z 2 ){i-\ zl \^^ 

7T 

Ml ^2 

(59) 

Further simplification of the obtained equations is possible if one assumes that 

1 < ^-^f^g^v.z). (60) 
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in the expansion Eq. (|5ip . This condition means a smallness of the source multiplied by the triple 
pomeron vertex function. Indeed, there is a suppression of this term compared to unity due to a 2 in 
front of Eq. (|59p if the source does not contain some large number (number of nucleons in nucleus, for 
example). In this case expanding the pomeron field Eq. (|57p one obtains 

$ M (y) = exp^y)^ - W\ X - 2lv fl B g (v,z , y) e 2 ^ (61) 

and the equation Eq. (i59j) becomes 

W\ l - 2lu fl B U^^y) e 2 ^u>, = J2 W^W^^)^^ (62) 

In Eq. ()62p we introduced a rapidity dependence in the function go as the price for this simplification 

g (v,Zo) -> g (v,z ,y) (63) 

since in Eq. (|62p this function possess some subdominant rapidity corrections. Thus in the case when 
Eq. (j62|) may be used instead Eq. (f59l) our ansatz becomes 



= exp(o;^)r MB 

kol 1 " 2 ^^ 50(^,^0,2/) (exp(^y)-l) + 1 1 ' 

The obtained equations are more complicated then that of Eq. (|44p in the high rapidity limit, because 
there is no factorization of z and v variables. 



6 The pomeron field at all rapidities and a transition region between 
small and large values of rapidity 

The obtained ansatzs for large and small values of final rapidity indicate, that the only function that 
will be changed due to the rapidity evolution will be a g function in the denominators of the ansatz. 
Therefore, as a definition of the large and small rapidities we could use a following observations. Between 
these two rapidity regions of the solution, at asymptotically large rapidity and small rapidity, exists a 
region of the transition between growing and saturated behavior of the pomeron field. In general, this 
region may be defined as a region where we can instead expansion in exponents e Wti y we use a asymptotic 
expansion in the exponents e _ajM y and vice versa. Therefore, from the forms of Eq. (|49p and Eq. (|57|) in 
this region may be defined by the following conditions 

\z cr \ l ~ 2%v f^ B (z cr ) g^ l (v) ( exp(w M y cr ) - 1 ) oc 1, (65) 

and 

WA 1 ' 21 " T^ B {z cr )go(v,z cr ) ( exp(w M y cr ) - 1 ) oc 1 (66) 

This, in analogy with the usual definition of saturation momenta, will define the critical scale of the z$ 
variable in the conformal basis 

(z cr ) oc g (v) (exp(u^y cr ) - 1) . (67) 

for which the transition occurs. It is interesting to note, that the source f^ B {z cr ) in Eq. (|67p brings 
some external scales dependence in Eq. (|67"|) . which defines the character of transition. We also could 
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consider the constraint Eq. (|67p as a definition of some "transition" rapidity y cr for fixed values of the 
final rapidity Y, vector zq and arbitrary external scales inside the source fn B 

y cr oc — In (l + g (u) \zcr\~ 1+2iv f^M ) . (68) 

Comparing Eq. (|65p and Eq. (|66p we see that one can define the transition region as a region, where the 
following relation is satisfied 

5oV) ~ 9o( v > z o)z =z cr (69) 

This relation shows, inter alia, that if we find the functional form of go(v, zq) , then this function must 
satisfy the convergent expansion in z$ around some value z cr determined by Eq. (|69p 

oo 

go(y, zq) = 9n X {v) ( z o ~ z cr) n ■ (70) 

n=0 

So, with the use of the g~o{v, Zq) function we could find the g^iy) function through Eq. (J7DJ) as well. 
Of course, for that we need to know a form of the go(u, zq) function, that for general case is not easy. 
Therefore, instead, we could define as a transition region the region where both terms of expansion 
Eq. (j6Tj) are equal 

exp(^y)f MB = \zv\ l - 2iv fl B Uv,^y) (71) 
Using the Eq. (|62|) we finally will obtain condition for the region where two solutions are overlapping 

e^f m (z cr ) = £ V^ 2 f, 1B ( Zl )f, 2B (z 2 )e^ + ^y-. (72) 

Clearly, this is a screening condition on the sources of the problem, which defines values of z cr and y cr 
for which the source of the projectile will be screened from the target by the triple pomeron interactions. 



7 The accuracy of ansatz 

The source of possible corrections to the pomeron field at large rapidities are the coefficients gk(y) for 
different k in the series expression for the pomeron field 

oo 

^(y,z) = z h - x z h ~ x 9k{v)e- k ^y (73) 

k=0 

In the similar expansion for the phenomenological pomeron the ration of the coefficients of the suc- 
cessive terms is proportional to a s . Based on this information we make an assumption that this is also 
the case in the expansion Eq. (I73p and we only need to find the overall normalization, i.e. the order in 
a s of the first term in the expansion. 

We plug this series into the equation of motion and write the first two equation from the chain of 
equations Eq. (I42p . In integrals over v we assume a contribution from such regions of v in which the 
BFKL structure of the series Eq. (|73p is kept, namely, from the "diffusion regions" of v with e Wv y ~ e^ y . 
In the first order of expansion in e~~ kuJl,y we obtain Eq. (|43|) 



w Ml gofa) = / di/2-| / dv 3 ^g (v 2 )g (v 3 )ni 2 (74) 



n 



One can easily see from counting the powers of a s that the coefficient function go (y) is of the order of 
l/a B . 

To show what are the limitations for the accuracy of the proposed ansatz we expand the expression 
in Eq. (J49]) for the large values of rapidity 

(/ \ / h-i-h-i \ 2 
1 + 1 g (v)\ +e- 2u »"{ 1 g (v)j + 

(75) 

or, approximately, 

z) = z h ~ l z h - 1 g Q {u) | 1 - e-«»y Z-Z^LZt go{v ) + e -^v ( ^f'' g Q (u)) + ... | (76) 

provided we neglect subleading terms in a s assuming 1 < 1 / a s . We see that the expansion Eq. (|76|) 
properly reproduces the expected (a^ 1 ) a; !') _1 behavior for each term, in the agreement with that of 
the phenomenological pomeron. 

Using this a s structure of the expansion Eq. ([TBI) one can easily see that the series converges for 
rapidities in the region of 

y > JL ln(J_) (77) 

a s a s 

Based on our discussion we conclude that the proposed ansatz Eq. (|64p has both the correct high energy 
behavior and the leading a s expansion in the rapidity region given by Eq. (I77D . If one wishes to take 
into account higher order corrections in a s in the expansion Eq. (|73p . the proper way to do this is to 
introduce a new form of the ansatz 

with some function 

oo 

F = E <fn{v,y) (79) 
n=-l 

where the all higher corrections are encoded. 

It is a straightforward procedure to show in a similar way that the same arguments also hold for the 
low rapidity expansion where one expands the ansatz in powers of e u ^ y . In this case the validity of the 
expansion is restricted to the rapidity region given by 

y < — H— ) (so) 

a s a s 

In order to account for higher order corrections in a s in the ansatz Eq. (|57|) one can also make a use of 
some function F in a way similar to that of Eq. (|78l) . However, in this case the function F possesses a 
different from Eq. (|79p form of expansion in the powers of a s 

oo 

F = £ o£ f n (u,y) (81) 
n=i 
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8 Conclusion 



In the present paper we discussed possible ways of an analytical solution to the BK equation in the 
conformal basis. We suggested the following ansatz for the solution of BK equation 



where the form of the function F in the denominator of Eq. (|82p depends on the region of rapidity where 
the solution is considered. The problem, therefore, is reduced to the evolution of this unknown function 
F with rapidity. Assuming the conformal invariance of the theory at high energy we simplify our problem 
proposing factorization of the coordinate dependence of the pomeron field in the this limit. This makes 
it possible to separate the coordinate and rapidity dependence resulting into Eq. ([301) . This equation is 
still not easy to solve, but it allows us to investigate the energy dependence of the solution to the BK 
equation not mixing it with the coordinate degrees of freedom. The next important simplification we 
make is keeping only the BFKL structure of the " fan" diagrams thus reducing Eq. (|3(jp to Eq. (|32|) . This 
stems from the well known fact that the BFKL propagators present the leading contribution at high 
energy in such diagrams. 

We find a solution to Eq. ([32]) which correctly describes high energy behavior of the exact solution 
and at the same time satisfies the initial condition at zero rapidity. Its expression is given by Eq. (|49p 
and has the same energy structure as the phenomenological "fan" amplitude. The matching between the 
correct high energy behavior and the fulfillment of the initial condition has also another aspect. It is 
related to the fact that a strong condition on the source given by Eq. ([54"]) is not satisfied in general, and 
thus dependence on a source breaks the conformal invariance of the solution. As one can easily see from 
the form of the solution Eq. ([49 1) at high energy the dependence on the source disappears restoring the 
conformal invariance. 

As a next step in our discussion we consider a question of a transition region in rapidity where a high 
energy solution ansatz transforms into the low energy one. This region can be thought of as one where the 
conformal invariance of the theory is restored or where the small rapidity expansion can be replaced by 
the asymptotic expansion. In this case some critical conformal scale can be introduced through Eq. (165[) - 
Eq. (j67p . Another way to find the behavior and the parametrical form of the critical scale is to match 
between of low and high energy ansatzs of the solution given by Eq. (|57p and Eq. ()49l) . respectively. In 
this case the critical scale z c r may be defined as a scale where the function ^(v, z) coincides with the 
function go -1 - 

It should be mentioned that the proposed ansatz is only an a approximation to the full solutions given 
by series Eq. (|4ip and Eq. ([5T]) Instead of using ansatz Eq. and Eq. (foTj) one can develop perturbative 
calculations and obtain a chain of equations similar to ones given by Eq. (I32p for large rapidities and by 
Eq. (|56p for of small rapidities. The calculation of coefficient functions of the expansion Eq. (|5ip valid for 
small rapidities corresponds to the calculation of the diagrams in the perturbative expansion. Therefore, 
the first equation ( Eq. (|56p or Eq. (|59p ) from a chain of equations for the coefficient functions is similar 
to ones obtained in Ref. [32J for the simplest "fan" diagrams. The only difference between our result 
and that of Ref. [32] is the conformal basis and normal coordinate representations correspondingly. The 
physical interpretation of the expansion Eq. (|41|) is not so clear. The expansion in the negative powers 
of exponents of rapidity cannot be put into one to one correspondence with diagrams. The high energy 
behavior of the solution to the BK equation in the coordinate representation is well known (see [35]) and 
the question of the relation between the conformal and coordinate representation at large rapidities will 
be addressed in our further studies. 

In this paper we considered the expansion of the coordinate dependent pomeron fields only in the 
conformal basis. Of course, for the practical applications the coordinate representation of the pomeron 




exp(^ M y) Tj 



(82) 
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field is more useful and convenient, but the expansion in conformal basis presents a more suitable frame- 
work for the investigation of the energy dependence properties of the pomeron field. These properties 
will determine the energy dependence of the full solution. The task of the inverse transformation of the 
found ansatz into the coordinate basis we leave , as we mentioned before, for further publications. 
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Appendix A: 



In this Appendix we calculate the integral which appears in Eq. (|23[) , To do this we want to rewrite it in 
terms of some dimensionless variables and use integral representation of the hypergeometric functions. 
First we rescale the variables Z\ and z 2 as follows 



oc 



Zfzf{z a - Z 1 )- A ^(Z 1 - Z 2 )-^{Z 2 - Z )~ A ™ dZ! dz 2 = (_l)-Ax2-A 25 ^i+72+2-A 6l -A 12 -A 2 



00 



-. N-A 12 -A n _7i+72+2-A 6l -Ai2-A 28 / 7l 72-A12 n (1 1 I /1 \-A ofi 7 , 

-1) 12 2 "z u u / w{w 2 (1 — wx) 01 ^1 J (1— W2) 20 dw\dw 2 



where u;i = £i/,Zo and ^2 = z 2 /zq. Now we perform the integration over 101. The relevant integral is 
given by 



h = r wT (1 - Wl y A ^ (l-—) 12 d Wl (A.2) 

./-oo V w 2 J 

As one can see this integral reminds the integral representation of the Gauss hypergeometric function 

C t b -\l - t)- b+c -\l - tz)- a dt = m l ( : C ~ h) 2F1 (a, b, c, z) (A.3) 
Jo T(c) 

but the care about the limits should be taken. We split the limits of the integration in Eq. (|A.3|) as 
follows 



00 rO pi poo 

= + / + / (A.4) 

-00 J —00 JO Jl 

The second term on the rhc in Eq. (|A.4[) is just the integral representation of the hypergeometric function 
given m Eq. (lA~3j) . The third term can be brought to the from of by substitution w = 1/t and in this 
case one obtains the following identity 

r t b ~\i - i)- 6+c - x (i - tz)- a dt = (-)-»- ^c-i r(o- c + i)r( c -&) 2F J c + 1 b + 1 i\ (A . 5) 

Ji r{a — o+l; V zj 

Thus the integral (0, 00) is obtained by summing Eq. (IA.3j) and Eq. (|A.5j) . The remaining part (— 00, 0) 
is obtained from the integral (0, 00) by substituting w = 1 — 1/t and reads 

° .b-ln ^-b+c-m ^-a^_/ x-c r(6)r(a-c) p / ^ „ n 1 



_ t )-^(i _ ^)-<^ = (_)-<= w v ' 2 iq a, c - 6, a - b + 1, (A.6) 

1 (a — c + + 1) \ 1 — 2/ 

Summing the expressions of Eq. ([A3]) . Eq. fA~5l) and Eq. (lA~6l) we can write the full expression to be 
used for w\ integration as follows. 



15 



£&-!(! _ t) -6+c-l (1 _ tz) -a dt = { _yc J-{b)T{a c) ^ & _ ^ _ ft + ^ + { ^ 

r(a-c + 6+l) V 1 — z, 

-a- 6+c -ir(a- c + i)r( c -6) p /_ _ , 1 _ , , , 1^ , r(6)r( c -6) 



+ ( _)-a- 0+ c-a^_^ / gF ao _ c+lo _ 6+1 + w v i 2 F X a,6,C,z) 

r(a — 6+1) V 2/ r(c) 

Comparing Eq. (|A.2|> and Eq. (|A.7|> we readily identify the parameters of Eq. (|A.7|) as 

o = Ai2, 6 = 1+71, c = 7i-A 5l + 2, z = — (A.8) 

Because of the inverse dependence of w 2 on z in Eq. (|A,8h it is more convenient for the further integration 
over w 2 to rewrite Eq. (|A.7|) in terms of the hyper geometric functions of the same argument -. This can 
be done using useful identities for the hypergeometric functions as follows ( see 15.3.4 and 15.3.7 in 

!) 



2 F 1 (a, 6, c; j^-) = (-l) a (l - zfz~ a 2 F 1 (a, c - b, c; 1) (A.9) 

and 

2 Fx(a,b,c;z) = (-l) a z- a l { ^lf - a \ 2Fl (a,l-c + a,l-b + a- -) (A.10) 

1 (6)1 (c — a) z 

^nrf"!! 2 Fi(M - c + 6, 1 - a + 6; i) 
r(a)T(c — 6) « 

With the help of Eq. (lA.9j) and Eq. (lA.lOj) the integral in Eq. (|A.7j) reads 

y°° t 6 - 1 (l-t)- 6+c " 1 (l-te)-°di = Ci- 2^1 6 -0+1,-0 + 6+1,^ + C 2 - 2-^1 (a,a-c+l,a-b + l, 
with the functions C\ and C 2 given by 

and 



zj r(c) 



. _ 6+c _ 1 r(a-c+i)r( c -6) , , lW /ivr( c -6)r(6- a ) , , 1V ^1 iy r(6)r(o-c) 



<* = r(a ; 6 Vl) +( - 1)a UJ r(c-a) ^-^l 1 r(a + 6-c+l) 

Thus with the help of Eq. ()A.8|) we identify the required integral Eq. (|A.2|) as 

Ii = Ci- 2i ? i(l + 7i,Ao 1 ,2 + Ai 2 + 7l ,u;2) + C 2 - 2 -Fi (a 12 , A i2 - 71 + A 5l - 1, A i2 - 71, (A. 11) 

with the functions C\ and C 2 given by 



d = {-l) l+ ^w. 



l+-y 1 „..l+7x r (l+7i)r(l-A 6l ) 

r( 7 i - A 6l + 2) 



2 
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and 



A,-, 1+ A 1 ,r(A 5l + A 12 - 7l -i)r(i-A 5l ) 1 , Al ^..A 12 r(i-A 5l )r(i + 7l -A 12 ) 



r(A 12 - 7l ) ' v ~ 2 r(2 + 72 -A 5l -A 



12 J 



+(-i ri -A 3l+2 (1 _ W2) ^ r(i + 7l )r(A 6l + A 12 - 7l -2) 

1 [A 5l + A i2 j 

The next step is to perform integration over variable w 2 . From Eq. (jA.ip with the definition of 
Eq. (IA.2j) we see that the integral over w 2 reads 



00 A 

hwf~ Al2 (1 - u; 2 )- A ^ du> 2 (A. 13) 

-oo 

It is clear from Eq. (|A.13|) and the result of the integration over w\ that the relevant integral is 

w a (l -wf 2 Fi (a,b,c,w)dw (A.14) 

As in the case of the integration over w\ we want to use the identity (see 7.152.5 in [35] ) 

( l w a (l - wf 2 F X (a, b, c, w) dw = ^g±j^f^ 3^2 (a, b, a + 1; c, a + /? + 2; 1) (A.15) 
Jo 1 (a + p + 2) 

and thus split the integration in Eq. (|A.14p as in Eq. (|A.4j) . The integral (1, oo) is obtained from Eq. (|A. 15|) 
by substitution w — > 1/w and with the help of the identity Eq. (jA.lOh . The result reads 

w a (l - wf 2 F X (a,b,c,w) dw = (A.16) 

-1 T + P— =r FTTvF? \ 3*2 (a, a - c + 1, a - a - (3 - 1; a - b + 1, a - a; 1 + 

r(a — a) r(o)T(c — a) 

, , 6+fl r(6-a-/9-l)r09 + l)r(c)r(a-6) „ „ , , „ , 

-1 b+/3 - ^ ^rTTFT TT3F 2 (b,b-c+l,b-a-P-l;-a + b + l,b-a;l) 

F(b — a) r(a)r(c — b) 

The sum of Eq. (1A.15|) and Eq. (|A.16h gives the contribution from the integration (0, oo). As in Eq. (|A.6j) 
the missing part (— oo, 0) is obtained from integral (0, 1) substituting w — > 1/(1 — w) and using identities 
Eq. CO]) and Eq. (lATTUl) 

o 

w a (l - wf 2 Fi (a,b,c,w) dw = (A.17) 

-oo 

T{a-a-/3- l)T(a + 1) T(c)T(b - c) 

-1 77 ^ FTTTFT V 3*2 (a, b - c, a - a - /3 - 1 -a - b + 1, a - 1) + 

T(a — pj r(o)T(c — a) 

. ^ n T(b-a- P-l)T(a + l)T(c)T(a-b) „ , , , „ „ , . 

(-1)^ - r(a)r(c-6) 3i?2 (-6 + c, a, 6 - a - /3 - 1; c, 6 - /9; 1) 

Finally summing the contributions of Eq. (|A.15p . Eq. (|A.16p and Eq. (|A.17p and plugging it into 
Eq. (|A.14p we obtain the contribution to Eq. (jA.ip coming from the integration over the holomorphic 
variables z\ and z 2 . The corresponding expression reads 



_ (-l) A oi~^B[l + 72 - A 12 , 1 - A 25 ]B[1 + A 12 , A 01 - A 12 ] (A.18) 
x 3^2 [{7i, " A oi' 1 + 72 - A 12 }, {-A12, 2 + 72 - A12 - A 25 }, 1] 



17 



_ (-i) A ei-Ai 2jB[1 + 72 _ Al2j 1 _ a 25 ]B(1 + Ai 2) A 5l - A 12 ] 
x 3^2 [{71, -Aqi, 1 + 72 - A 12 }, {-Ai 2 , 2 + 72 - A12 - A 25 }, 1] 

_ (-i)- A o 1 +Ai 2jB(1 + 72 _ Al2 , 1 - A 25 ]5(l + Ai 2 , A 8l - A 12 ] 
x 3 i ? 2 ({7i, -Aqi, 1 + 72 - A12}, {-A12, 2 + 72 - A12 - A 20 }, 1] 

_ (_i)72-A 5lB(1 + 72 _ Ai2> _i _ 72 + Al2 + A 25 ]£(l + Ai 2 , A 6l - A 12 ] 
x 3^2(171, -A 5l , 1 + 72 - A12}, {-A12, 2 + 72 - Ai 2 - A 25 }, 1] 

_ (_ 1) 72+A 0l -2A 12jB(1 _ 7i + 72 _ Ai2> _! + 7l _ 72 + Al2 + A 2 - ]B(1 + 7l + Ai 2 , A 5l - A12] 
x 3^2({7i, 1 + 7i + A12, -1 + 71 - 72 + A12 + A 20 }, {1 + 71 + A 5l ,7i - 72 + A 12 }, 1] 

(_l)7i-A ( - n -2( 72 -A 12 )-A 12 -A 20j g ( _ 7i _ A6i) ! + 7l + Al2 ]^(l + A12, -Ai 2 ]fl(-1 + 71 - 72 + A12 + A 25 , 1 - A 26 ] 

B(l - 71,71] 

x 3^2({7i, 1 + 7i + A12, -1 + 71 - 72 + A12 + A 2 g}, {1 + 71 + A 5l ,7i - 72 + A 12 }, 1] 

(_l)7i-A ( -, 1 -2( 72 -A 12 )+A 12 -A 2( - lj g ( _ 7i _ A . i; Agi _ Al2 ]g(l + A12, -A 12 ]B(-1 + 71 - 72 + A 12 + A 25 , 1 - A 26 ] 

£(1 + A 5l ,-A 5l ] " ' 
x 3 ^2({7i, 1 + 71 + A12, -1 + 71 - 72 + A12 + A 25 }, {1 + 71 + A 6l , 7 i - 72 + Ai 2 }, 1] 

_ ( _ 1) - 71+ A 5l -A 12 -A 25B(1 + ^ + Ai2) _ Al2 ] B (_! + 7l _ 72 + Al2 + A ^ 1 _ A25 ] 

x 3 *2({7i, 1 + 7i + Ai 2) -1 + 71 - 72 + A12 + A 25 }, {1 + 71 + A 5l , 71 - 72 + Ai 2 }, 1] 

_ (-l)72+A 5l -2A 12£ > (A5i _ Ai2) 2 + 72]jB (_l + 7l _ 72 + Al2> 1 + 72 _ Al2 ] 

x 3 F 2 ({2 + 72, 1 + 72 - A12, A 20 }, {2 - 71 + 72 - A i2 , 2 + 72 + A 5l - A 12 }, 1] 

^ (-1)-A 5lj g(2 + 72, 1 - A 26 ]fl(l + A12, -A12] 



(l + 7i + Ai 2 )B(7i,2 + Ai 2 ] 
x 3 F 2 {{2 + 72, 1 + 71 + A12, 1 - A 6l + A12}, {2 + A12, 3 + 72 - A 20 }, 1] 



(_l)72-A 5lB(2 + 72) _ 2 _ 72 + A 25 ]B(1 + A12, -A 
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(l + 7i + Ai 2 ) J B( 7 i,2 + A 12 ] 
x 3 F 2 ({2 + 72, 1 + 71 + A12, 1 - A 6l + A12}, {2 + Ai 2 , 3 + 72 - A 25 }, 1] 

+ (-l) A oi^(2 + 72, 1 - A 26 ]B(1 + 71 + A12, -1 - A12] 
x 3 F 2 ({2 + 72, 1 + 71 + A12, 1 - A 6l + A12}, {2 + A12, 3 + 72 - A 20 }, 1] 

( _ 1) -2A 6i -2( 72 -a 12 )-a 12 -a 26 B{1 + Ai2) - Al2 ]B(-l - 72 - A 6l + A 12 + A 26 , 1 - A 25 ] 

(7i + A 6l )B( 7 i,l + A 5l ] 
x 3 F 2 ({-A 5l , 1 - A 6l + A12, -1 - 72 - A 5l + A12 + A 20 }, {1 - 71 - A 5l , -72 - A 6l + A 12 }, 1] 
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( _ 1) -2A 6l -2( 72 -A 12 )+A 12 -A 26j g (1 + Al2> _ Al2 ]j?(-l - 72 - A 6l + Aig + A 2 g, 1 - A 2g ] 

( 7l + A 6l )B( 7l ,l + A 5l ] 
x 3 F 2 ({-A 5l , 1 - A 6l + Ai 2 , -1 - 72 - A 6l + A12 + A 20 }, {1 - 7l - A 5l , -72 - A 6l + A 12 }, 1] 

- (-l) 72 - Al2 S(A 5l - A12, -1 - 72 - A 5l + Ai 2 + A 26 ]B(-1 + 7l - 72 + A12 + A 25 , 1 + 72 - A 12 - A 20 ] 
x 3 F 2 {{A 2 - , -1 + 7l - 72 + A12 + A 25 , -1 - 72 - A 6l + A12 + A 25 }, {-1 - 72 + A 25 , - 72 + A 12 + A 25 }, 1] 

(_1)72-a (1i -a 12B(1 + Ai2> _ Al2 ]ff(-i + 7l _ 72 + Al2 + a 26; 1 + 72 _ Al2 _ a 25 ] 

(-1 - 72 - A 5l + Ai 2 + A 25 )B(1 - A 5l + Ai 2 , -1 - 72 + A 25 ] 
x 3 F 2 ({A 2 - , -1 + 7l - 72 + A12 + A 25 , -1 - 72 - A 01 + A12 + A 25 }, {-1 - 72 + A 2fj , - 72 + A 12 + A 25 }, 1] 

_ (-lp-Am^x + Al2 , -A 12 ]ff(-1 - 72 - A 6l + Aig + A 25 , 2 + 72 - A 20 ] 
(-1 + 7i - 72 + A12 + A 25 )B(7i, -72 + A12 + A 25 ] 
x 3 i ? 2 ({A 25 , -1 + 71 - 72 + Ai 2 + A 25 , -1 - 72 - A 6l + A 12 + A 25 }, {-1 - 72 + A 25 , - 72 + A 12 + A 25 }, 1] 

The integration over antiholomorphic variables z\ and z 2 now can be easily performed using this last 
result. 
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